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We study the dynamics of the front separating a spatio-temporally chaotic region from a stable 
steady region using a simple model applicable to periodically forced systems. In particular, we 
investigate both the coarsening of the front induced by the inherent 'noise' of the chaotic region, 
and the long wavelength dynamics causing the front to develop cusps. 



INTRODUCTION 

In this paper we study the dynamics of the front sep- 
arating a spatio-temporally chaotic state from a stable 
steady ordered state. Such situations occur in many ex- 
perimental settings. In an experiment on a vertically 
vibrating granular monolayer of spheres [l[ both a state 
at rest on the plate and a chaotically bouncing state are 
stable. When a small perturbation is applied to the sta- 
tionary state, the chaotic state is observed to invade the 
stable state through a propagating front. In a Rayleigh- 
Benard convection experiment Q both straight rolls and 
spiral defect chaos are stable under some conditions and 
it is observed that a region of straight rolls is invaded by 
a region of spiral defect turbulence. 

For our study we employ a type of model called a Con- 
tinuum Coupled Map (CCM) introduced in Ref. 
Models of this type (Sec. II) are appropriate to peri- 
odically forced systems (such as that in the experiment 
of the Ref. (H). In common with other generic mod- 
els, like the complex Ginzburg-Landau equation or the 
Swift-Hohenberg equation, CCM models are meant to 
incorporate the minimal basic properties capable of re- 
producing the phenomena of interest. With this in mind 
we construct our CCM model to incorporate the essen- 
tial feature that both a stable steady homogeneous at- 
tractor and a spatio-temporally chaotic attractor exist. 
Using our CCM model, we numerically investigate two 
phenomena: (i) the coarsening of the front due to the in- 
herent 'noise' associated with the spatio-temporal chaos 
(Sec. Ill), and (ii) cusp formation induced by initial long 
wavelength perturbations of the front location from the 
fiat state (Sec. IV). 

With respect to (i), an important concept used to 
study various coarsening processes is scaling. For a large 
number of systems (e.g., see Ref. Q), it is found that 
the interface width due to roughening, w(t), increases as 
a power of time, w(t) ~ t 13 . The width eventually satu- 
rates at a value that increases as a power of the system 
size, w{L x ) ~ L". These scaling properties are also ob- 
served in our model, and we determine and discuss the 
scaling exponents a and /3 that we find. 

With respect to (ii), we argue that on long scale (i.e., 



long compared to w(t)), our fronts propagate at constant 
speed in a direction locally normal to the interface. We 
show that this basic property explains the mechanism 
of cusp formation and the evolution of the shape of the 
fronts observed in our numerical simulations. 



MODEL 

As in reference we consider a CCM model which 
maps a field £ n (x) forward from time n to time n + 1. 
With reference to a system driven periodically in time 
(e.g., as in [l[), we may think of £ n (a?), with n = 1, 2, • • • , 
as being the system state stroboscopically sampled once 
each period. Furthermore, we consider x = (x, y) to be 
two-dimensional, and, for simplicity, we take to be a 
scalar field. The CCM model mapping £ ra to consists 
of two steps: The first step is a nonlinear local operation 
in which a one-dimensional map M is applied to (a?) at 
each point in space, 



C(x) = M(Ux))- 



(1) 




FIG. 1: Schematic (not to scale) of the model map. 



2 




FIG. 2: Spatio-temporal chaos of CCM model for A = 7.0, r = 0.4, and w = 0.29: (a) A snapshot of fji(x). Bright regions 
indicate large amplitude, (b) Wavenumber power spectrum, averaged over 100 frames. £(fc) is the Fourier transformation of 
£(a?). (c) Time correlation function C(r), Eq. ([7]). 



The second step is a translationally invariant linear oper- 
ation coupling the dynamics at nearby spatial locations. 
The most general such coupling is conveniently expressed 
in terms of the spatial Fourier transform. If is the 

spatial Fourier transform of £„(:?), then we write £ n +i(fc) 
as 



L + l(k) = f(k)i' n (k), 



(2) 



from which £ n+ i(x) is obtained by inverse Fourier trans- 
forming. 

The model is then specified by the choice of the non- 
linear map M and the linear spatial coupling f(k). We 
make these choices so as to include the minimum prop- 
erties that we hypothesize are relevant for the investi- 
gated phenomena Since we desire the simultaneous 
existence of a stable steady state as well as a spatio- 
tcmporally chaotic state, we choose the map M to have 
a stable fixed point attractor and a chaotic attractor. A 
convenient choice having this property is given by 

M(0 = r£ + Aexp[-(£ - \flo 2 \ - A=xp(-l/<7 2 ). (3) 

Referring to Fig. HJ we see that this map has a stable 
fixed point at £ = 0, for r < 1. Moreover, any initial 
point in £ < u is attracted to this point. We also see from 
Fig. [T]that the interval u < £ < v is mapped into itself. 
Thus there is (at least) one attractor in this interval. 
For the parameter values A, r, and w that we investigate 
(^4 = 7.0, r = 0.4, (7 = 0.29) there is one attractor in 
u < £ < v and it is chaotic. 

Our choice of the spatial coupling f(k) is similarly mo- 
tivated by a desire for simplicity. We assume that the 
coupling is isotropic. Thus we can write / as f{k) = /(fc), 
where k = Taking f(k) > we write 



where "f(k) is a wavenumber-dependent growth/damping 
rate per period. Since we want the spatio-temporal chaos 
to have a finite spatial correlation scale, such a scale must 
Ik- reflected in our choice of f(k). Denoting this scale by 
k$ , we make the following simple choice |3[ for 'y(k), 



7 (fc) = 



2 V fc, 



1 - 



(5) 



/(*) = exp[ 7 (fc)], 



(4) 



Thus 7(fc) > (growth) for fc < fco, l{k) has its maxi- 
mum value at fc = fco, and "f(k) becomes strongly negative 
(damping) as k becomes large. 

Our numerical implementation of this CCM model cm- 
ploys doubly periodic boundary conditions with period- 
icity lengths L x in x and L y in y. The nonlinear map 
operator is applied at points on a square grid, while the 
spatial coupling operator ^ employs fast Fourier trans- 
forming from x to k and back. 

That £(a?) = is an attractor can be seen by introduc- 
ing an initial perturbation at wavenumber fc, i5£ocxp(ifc ■ 
x). Linearization of the CCM model about £o(af) = 
then shows that this perturbation evolves with time to 
5£„exp(i£ • £), where 8£ n = 5£ [M'(0)f(k)} n ,M'(0 = 
dM/dS,. For the parameters we choose M'(0)/(fc) < 1 
for all k (in particular M'(0)f(k Q ) < 1). Thus the ho- 
mogeneous state £(:?) = is an attractor for the system. 
We also find that, as we had anticipated, for other ini- 
tial conditions there is another attractor which is spatio- 
temporally chaotic. 

Figure [2] shows the properties of the spatio-temporal 
chaos produced by our model. FigureElJa) shows the spa- 
tial pattern £, n (x) at a representative time. This picture 
applies to a time n = 45 evolved from an initial con- 
dition where £o(a?) was chosen randomly with uniform 
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distribution between £o = and £o = 7.5. Visually, we 
observe that the pattern appears to have a characteris- 
tic scale of the order of k^ 1 . This is confirmed by the 
wavenumber power spectrum, Fig. [2{b). We note that 
the only length scales in our model are k^ 1 , the system 
size L ~ L x ~ L y , and the grid size 5, and that, by our 
choice L ^> k^ 1 » 8, we had sought to obtain spatio- 
temporal chaos with properties independent of L and 6. 
Figure EJb), which evidences variation on the scale fco, 
conforms with this expectation. Further discussion of 
the form observed for |£(/c)| 2 is given in the appendix. 
To characterize the temporal variation of the patterns, 
Fig. [2jc) shows a plot of the time correlation function 
C(t) defined as 



C(r) 



1 



N~N, 



1,3 



<es)-zh> 



(6) 



, (7) 



where < ... > means time average, £ =< £ >,N X}V = 
L x y /S, and denotes the (x,y) location of a grid 

point. As can be seen from Fig. &[c) the time correlation 
function decays to zero with increasing time r (where r 
is an integer), confirming that the temporal behavior is 
chaotic. 

We have also examined other parameter values for 
which CO)-© yields spatio-temporal chaos, and we find 
behavior similar to that in Figs. [SJ 



PROPAGATION OF A FLAT FRONT 

The main objective of our investigation is to charac- 
terize the dynamics of the interface between a spatio- 
tcmporally chaotic state and a stable steady state. In 
our first set of simulations we focus on an initially flat 
interface, y = yo at t = 0. After we generate an initial 
spatio-temporally chaotic state (Fig. E^a)), we create the 



y-y n 
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interface by setting the amplitude £ of all grid points with 
y < yo to the stable steady state value £ = 0. During fur- 
ther iterations, the front between the chaotic and steady 
state moves downward, i.e., the chaotic state propagates 
into the stable steady state, and the front coarsens (see 
Fig. [3]). In order to examine the front dynamics for 
long times and save computational time, we use a shift- 
ing method: On every iterate, we reset £ to zero in the 
region adjacent to the bottom (y = 0) of the periodic 
box, and, when the front comes close to y = 0, we shift 
the whole system upward in the y direction. Due to the 
periodic boundary conditions, after this shift there will 
be a region below the front and above y — that is in the 
spatio-temporally chaotic regime, and we then set £ = 
in this region. 

After an initial transient, the scaling properties of the 
front are studied using the following definition of the in- 
terface width. First, we calculated the average value of £ 
at fixed y, £(?/) = J Q Lm £(x,y)dx/L x . We then note that 
the basin boundary between the two attractors of the 
one dimensional map M is at the unstable fixed point, 
u = 0.4824 (Fig. [TJ), and that the average £ for the 
spatio-temporally chaotic state (Fig. Ufa)) is approxi- 
mately four times the critical value. Thus, we define a 
lower boundary of the front, yi, by = u and an up- 
per boundary of the front, ?/2, by £(2/2) = 3it. The width 
of the front is then defined as (2/2 — Ui)- 

Because of the inherent 'noise' generated in the spatio- 
temporally chaotic region, the proper quantity to study 
is the ensemble averaged mean of the front widths. We 
calculate ensembles using many different random initial 
conditions. Our results for the ensemble averaged width 
w of the front arc obtained by averaging (j/2 — y\) over 
10 runs for the largest system L x = 1024 and over 500 
runs for the smallest L x = 64. The typical coarsening of 
the front is shown in Fig. [3l 

One observation from our simulations is that the prop- 
agation velocity of fronts is constant, except for a few 
transient initial iterations. That is, the velocity does not 
depend on time or system size. A constant propagation 
front velocity is also observed in the experiments in Refs. 
i and 0. To minimize the effect of the initial transient, 
we redefine time as the total increase in the area of the 
chaotic state. 

As is typical for a front coarsening problem the 
time and system size dependence of the mean front width 
w can be described by a scaling function g{u), 



(8) 



FIG. 3: Coarsening of a flat front. L x = 1024. 



Here g(u) is constant for u <C 1 and g{u) as vT^ for 
u ^> 1 . For t 3> L x , the width saturates at w ~ L x where 
a — [3z is the roughness exponent. Barabasi and Stanley 
[1] have summarized the values of the scaling exponents 
z, a, and (3 that are obtained for several experimental 
systems as well as relevant theoretical results. 
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In Fig. BJa). we show w versus t for different system 
sizes. These data show two characteristic regimes: power 
law growth, followed by saturation. The growth expo- 
nent, /3, is calculated by measuring the slopes of straight 
line fits to the data before saturation, and the roughness 
exponents a is calculated by comparing the saturation 
widths. Figure Ub) shows w/L" versus t/L*, where a 
and z have been adjusted to a = 0.49 and z = 1.81 (cor- 
responding to /3 = a/z = 0.27). Consistent with Eq. 
((HI), we find collapse of the data in Fig. IUa) to a sin- 
gle scaling function. The exponent values we obtain are 
roughly consistent with those of both the two dimensional 
Eden model (a = 0.5, /3 = 0.3) and the two dimensional 
Kardar-Parisi-Zhang equation {a = 1/2, /3 — 1/3) 0, Sj- 
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FIG. 4: Scaling of width with the size: (a) unsealed, (b) 
scaled with a = 0.49, z = 1.81 from which /3 = a/z = 0.27. 
Circles (L x — 64), diamonds (L x — 128), squares (L x = 256), 
triangles (L x — 512), and stars (L x — 1024). L y — 256 for all 
cases. 
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FIG. 5: Geometrical picture of the trajectories given by Eqs. 
(fTO)) and |TT]) (a) for t < t c and (b) for t > t c . 

EVOLUTION OF A NON-FLAT FRONT 

We now consider the evolution of a front on large 
length scale. Specifically, we are interested in the case 
where the front is initially not flat; that is, the position 
of the front is initially given by yo = h(xo). Further- 
more, we assume that, as the front evolves, the scale on 
which the front position varies, I w (h/h'), remains large 
compared to the front width, 



/ > w. 



(9) 



To analyze this situation, we consider that a point on 
the front moves with a normal velocity v whose magni- 
tude, \v\ = v, is constant in time. This assumption also 
implies that the direction of v following the trajectory of 
a point on the front is constant in time. This is because 
the slope of the front, dx/dy, following a trajectory does 
not change with time; i.e., it depends only on the initial 
location xq on the front and not on t. This is illustrated 
by the construction in Fig. Eta). As shown in the figure, 
the trajectory line segments are straight, are all of the 
same length, vt, and are normal to both the initial front 
and to the evolved front. Considering an initial front po- 
sition given by yo = h(xo), propagation at a velocity v 
normal to the front then yields 



x{t) 

y(t) 



x -{vt)h'(x )[l + {h'(x )) 2 } 
h(x )-vt[l + (h'(x )) 2 )- 1 / 2 . 



-1/2 



(10) 
(11) 
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FIG. 6: Comparison between theoretical curves and our 2D 
simulations for two different system sizes. The front curve 
from our model is defined as the smallest y value (for given 
x) at which £ = u, where u = 0.4824 is the basin boundary 
point depicted in Fig. [T] 

At any given time Eqs. (|10|) an d (fTTj) specify the front 
position parametrically with xq as a parameter. As an 
example consider the case of an initial sinusoidal undu- 
lation of the front yo = h(xo) = C cos(kxo) . As the front 
propagates, the initial sinusoidal curve becomes distorted 
so that the maxima become sharp and the minima be- 
come broad. As can be seen from Fig. [5] (a) this arises 
because of the converging (diverging) of trajectories that 
originate near maxima (minima). 

Past a critical time t = t c a cusp develops at the max- 
ima of the evolved front Q . The cusp formation time is 
determined by noting that dx/dxo first becomes zero at 
t = t c . From Eq. (JTUJ) we obtain 



For t > t c there are pairs of values of Xq for which the 
trajectories given by (fl"0|) and (fTTj) pass through each 
other. For a given time t greater than t c we refer to the 
range of xq for which this occurs as the unphysical range. 
The development of the unphysical range is illustrated in 
Fig. [5] (b) where the dashed portion of the curve shows 
the result of plotting (fT0|) and (jTTJ) for the unphysical 
range. 

Figure [5] shows how cusps develop in time. In Fig. [5] 
the noisy front curves are from our simulations and the 
smooth curves are from Eqs. (fT0|) and (fTT|) with x re- 
stricted to the physical range. The two cases shown in 
Fig. ffl (namely L x = L y = 128 and L x = L y = 1024) 
illustrate how front roughening becomes of less influence 
as (HJ becomes better satisfied. For both cases in Fig. [6] 
the initial sinusoid has amplitude kC = 1 and wavenum- 
ber k = 4ir/L x . Thus from (fl"2|) we have t c ~ 1/k ~ L x . 
At times t ~ t c , the roughening is a small effect if kw(t c ) 
is small. Since w ~ (assuming t/L* < 1), we see that 
kw(t c ) ~ l/L^T 13 , and roughening will be inconsequen- 



tial for the large scale front evolution if L x is sufficiently 
large. The good agreement of the L x = 1024 numerical 
results from our CCM model with the theory, Eqs. (fT0|) 
and (fTTj) . confirms that the front does indeed propagate 
at constant velocity in a direction perpendicular to the 
interface. In particular, for long length scales I ^> w 
examined in Fig. [6] we see no evidence for curvature de- 
pendence of the front velocity. 
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FIG. 7: The map Eq. JT3J. 



As a comparison, we have also considered the propaga- 
tion of a front between two steady homogeneous states. 
In particular, replacing M(£) in Eq. [3] by 

M(f) = 0.1 + tanh(2£), (13) 

we see (Fig. [7]) that there are two stable fixed points, 
a and /3, and one unstable fixed point u. Using (|13|). 
([2"| and(0 with a sinusoidal front with initialization of f3 
above the front and at a below the front, we find that 
the front shape evolves smoothly (with no roughening) 
according to Eqs. fTU|) and (fTTj) . Thus the cases of a 
chaotic invading region and a nonchaotic invading region 
become similar for large L x (e.g., L x » w » 1/feo)- 

In conclusion, we have introduced a continuum coupled 
map model for the study of the dynamics of a front sep- 
arating a region of spatio-temporal chaos from a stable 
steady region. This model is applicable to periodically 
forced systems. We find that the front roughens and 
that this coarsening obeys a scaling hypothesis, Eq. ©. 
We also investigate the large length scale evolution of a 
nonplanar front. We find that this evolution is consistent 
with the hypothesis that, on large scale, the front veloc- 
ity is constant and normal to the front. This hypothesis 
and our numerical simulations indicate the formation of 
cusp structures in the front. 

This work was supported by ONR (Physics). The 
work of SCV was supported by a research fellowship from 
the Alfred P. Sloan Jr. Foundation and by the NSF 
(DMR9975533). 
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which is plotted in Fig. [8] as the clashed line along with 
the data from Fig. [2](b) . It is seen that Eq. (A.l) pro- 
vides a crude indication of the general form of |£(fc)| 2 . 
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FIG. 8: |£| versus k compared with Eq. (A.l). 



APPENDIX 

We now comment on the specific form that we have 
found for |£| 2 (Fig. [2£b)). In this connection we note 
that in the limit of a wildly varying map M(£) with Lya- 
punov exponent approaching infinity £J n (x) will be wildly 
varying in space. This is because small variation of £, n {%) 
with x are greatly amplified when M is applied. Thus, 
in this limit, the spatial correlation function for £' n {x) 
will be a delta function, and |£'„(fc)| 2 =< (€' n (x)) 2 > 



independent of k. Thus from {2| 

m\ 2 =< (cm 2 > f(k), 
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